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ABSTRACT 

We present a dynamical study of the Galactic black hole binary system Nova Muscae 1991 (GS/GRS 
1124—683). We utilize 72 high resolution Magellan Echellette (MagE) spectra and 72 strictly simultane¬ 
ous V-band photometric observations; the simultaneity is a unique and crucial feature of this dynamical study. 

The data were taken on two consecutive nights and cover the full 10.4-hour orbital cycle. The radial veloci¬ 
ties of the secondary star are determined by cross-correlating the object spectra with the best-match template 
spectrum obtained using the same instrument configuration. Based on our independent analysis of five or¬ 
ders of the echellette spectrum, the semi-amplitude of the radial velocity of the secondary is measured to be 
K 2 = 406.8 ±2.7 km s^', which is consistent with previous work, while the uncertainty is reduced by a factor 
of3. The corresponding mass function is/(M) = 3.O2±O.O6M0. We have also obtained an accurate measure¬ 
ment of the rotational broadening of the stellar absorption lines (vsini = 85.0 ±2.6 km s^*) and hence the mass 
ratio of the system ^ = 0.079 ± 0.007. Finally, we have measured the spectrum of the non-stellar component of 
emission that veils the spectrum of the secondary. In a future paper, we will use our veiling-corrected spectrum 
of the secondary and accurate values of K 2 and q to model multi-color light curves and determine the systemic 
inclination and the mass of the black hole. 

Subject headings: black hole physics — stars; black holes — binaries: general — X-rays: binaries 


1. INTRODUCTION 

Mass is the fundamental parameter for a black hole. Ac¬ 
cording to the No-Hair Theorem (e.g., Israel 1967; Hawking 
1971), mass and spin together provide a complete description 
of an astrophysical black hole. An accurate measurement of 
mass is a prerequisite for measuring black hole spin via the 
continuum-fitting method (see McClintock et al. 2014 for a 
review). Stellar-mass black holes in the Milky Way and neigh¬ 
boring galaxies are discovered in X-ray binary systems, some 
of which are persistent X-ray sources, while others are X-ray 
transients which have gone through one or more X-ray out¬ 
bursts (Remillard & McClintock 2006). About two dozen of 
these black holes have been dynamically confirmed to have 
masses in the range M = 5-30 Mq. 

X-ray Nova Muscae 1991 (hereafter NovaMus), with an or¬ 
bital period of 10.4 hr, is one of about a dozen black-hole tran- 
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sient systems that are distinguished by their short orbital peri¬ 
ods, P < 12 hr. (For a sketch to scale of nine of these systems, 
including the prototype A0620--00, see Figure 1 in McClin¬ 
tock et al. 2014.) NovaMus was discovered during its 1991 
outburst independently by the Ginga (Makino et al. 1991) 
and GRANAT (Lund & Brandt 1991) missions. Eight days af¬ 
ter discovery, it reached a peak X-ray intensity of 8 Crab in 
the 1-6 keV band (Ebisawa et al. 1994). Over the course of 
eight months, the source passed through all the canonical X- 
ray states before returning to its quiescent state (Ebisawa et al. 
1994; Remillard & McClintock 2006). Esin et al. (1997) used 
the 1991 outburst data for NovaMus to develop a model of the 
states that combines the standard model of a thin accretion 
disk (Shakura & Sunyaev 1973) and the advection-dominated 
accretion flow (ADAE) model (Narayan & Yi 1994). 

Interestingly, hard X-ray observations of NovaMus made 
at the peak of the outburst revealed evidence for positron- 
electron annihilation in the form of a relatively narrow and 
variable emission line near 500 keV (Sunyaev et al. 1992; 
Goldwurm et al. 1992). Kaiser & Hannikainen (2002) argued 
that the line could be generated by annihilation in the bipo¬ 
lar outflow of the system. Martin et al. (1996) suggested an 
alternative mechanism, namely, Li production near the black 
hole that gives rise to a 476 keV nuclear line; their support¬ 
ing evidence is their detection of a relatively strong lithium 
absorption line A6708 in the optical spectrum. 

Soon after its discovery, NovaMus was considered to be a 
strong candidate for a black hole binary based on its X-ray 
spectral properties and light curve, and its similarities to the 
prototype black-hole X-ray transient A0620—00 (McClintock 
8c Remillard 1986). Eollowing the discovery of the optical 
counterpart (Della Valle et al. 1991), the black hole nature of 
the primary was established in the customary way by measur¬ 
ing the mass function: 
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where K 2 is the semi-amplitude of the velocity curve of the 
secondary, i is the orbital inclination angle of the binary and 
q is the ratio of the companion star mass M 2 to that of the 
compact primary M. Three dynamical studies performed in 
quiescence found a value for the mass function (the minimum 
mass of the compact primary) of « SMq (Remillard et al. 
1992; Orosz et al. 1996; Casares et al. 1997), which is widely 
considered to be the maximum stable mass of a neutron star 
(Kalogera & Baym 1996). 

Orosz et al. (1994) obtained a quite uncertain estimate of 
the mass ratio q = M 2 /M by measuring the radial velocity of 
the Ha emission line. Casares et al. (1997) made a conven¬ 
tional and direct measurement of the mass ratio by measuring 
the rotational broadening of the photospheric absorption lines 
of the secondary; they obtained a result that is consistent with 
that of Orosz et al. (1994). 

The inclination angle i is usually estimated by modeling 
the ellipsoidal modulation of the light curves of the secondary 
star after a system has returned to its quiescent state. This 
modeling is complicated in the case of the short-period sys¬ 
tems because the stellar component of light is significantly 
contaminated by the emission from the accretion disk; this 
so-called “disk veiling,” can be quite significant (> 50% in 
the T-band) and variable in “active” quiescent states (Cantrell 
et al. 2010). Orosz et al. (1996) estimated a disk contribution 
at ~ 5000 A that could be as high as 50%. Meantime, using 
different data, Casares et al. (1997) estimated a disk contri¬ 
bution at Ha A6563 of 15%, suggesting that disk veiling may 
be less significant at longer wavelengths. Based on this no¬ 
tion, Gelino et al. (2001) modeled near-infrared {J and K) 
light curves of NovaMus assuming that the disk contribution 
could be ignored in these bands, and they thereby estimated 
the systemic inclination to be i = 54° ± 1.5°; their estimate of 
black hole mass (which also depends on the mass ratio q) is 
M = 6.95 ± 0.6 Mq. 

These estimates of i and M are quite uncertain because they 
rest on the assumption that the effects of veiling are negligible 
in the near-infrared, while there is considerable evidence that 
it may not be (e.g., Reynolds et al. 2007, 2008; Kreidberg et 
al. 2012). In this paper, we obtain the first detailed and secure 
measurement of the spectrum of the non-stellar component of 
light that veils the stellar spectrum. While this achievement is 
in part due to the quality and quantity of our data, it is based 
primarily on the strict simultaneity of our spectroscopic and 
photometric observations. The simultaneity is crucial because 
of the aperiodic variability of the non-stellar component. Our 
two other principal results are a greatly refined measurement 
of the /T-velocity and the first robust measurement of the mass 
ratio q, which differs significantly from the earlier estimates 
mentioned above. 

This paper is structured as follows. Observations and data 
reduction are described in ij2] We derive the mass function via 
a radial velocity analysis in S|3] The mass ratio is determined 
via rotational broadening in S|4] The disk veiling is also mea¬ 
sured in this section. The light curve of NovaMus is presented 
in ^ Finally, we discuss our results in ij 6 ] The light curve 
modeling and the determination of the inclination and black 
hole mass will be presented in a subsequent paper. 

2 . OBSERVATION AND DATA REDUCTION 
2.1. Magellan Echellette Spectroscopy 

The spectroscopic data were obtained using the Magel¬ 
lan Echellette spectrograph (MagE; Marshall et al. 2008) 


mounted on the 6.5-meter Magellan/Clay Telescope at the Las 
Campanas Observatory (ECO) in Chile. All the data were 
taken on the nights of 2009 April 25-26 UT. A total of 72 
spectra were obtained, each with an integration time of 600 s. 
During the same two nights, and with the same setup, we ob¬ 
tained spectra of 38 stars, principally K dwarfs, that serve 
as radial velocity templates. The seeing was similar on both 
nights, with a typical value of sa 0.7" and a range of 0.5"-l .3". 
Each MagE spectrum consists of 15 orders (corresponding to 
sequence numbers #6-20). The full wavelength coverage of 
the instrument is 3100-11000 A, with order #20 at the blue 
end (central wavelength 3125 A) and order #6 at the red end 
(central wavelength 9700 A). The spatial scale is « 0.3"/pixel. 
The spectral dispersion increases from 0.2 A per pixel for or¬ 
der #20 to 0.7 A per pixel for order # 6 , while the velocity dis¬ 
persion is 22 km s^* per pixel for all orders. For our chosen 
0.85" slit, the spectral resolution is « 1 A at 5000 A, corre¬ 
sponding to « 60 km s^*. 

The MagE data were reduced using the pipeline developed 
by Carnegie Observatorie0> which performs bias correc¬ 
tion, flat-fielding, and wavelength calibration. The bias was 
corrected using the overscan region of the detectors. Three 
sets of flats were utilized in the calibration: 1 ) order defini¬ 
tion flats taken with the Xe-Flash lamps; 2) blue flats to de¬ 
fine the flat field at the blue end of a spectrum (also taken 
with the Xe-Flash lamps); and 3) red flats to define the flat 
field at the red end of a spectrum and to correct the fringing, 
which were taken with the Qf lamps. The wavelength calibra¬ 
tion was performed using arc spectra obtained with Thorium- 
Argon lamps. The RMS residual of the wavelength solution 
is < 0.05 A for orders #8-18 and < 0.1 A for the remaining 
orders. The performance of wavelength calibration was exam¬ 
ined using sky emission lines; the relative wavelength offset 
is negligible. 

The mean MagE spectrum of NovaMus in the observer’s 
rest-frame is shown in Fig.[T] which covers spectra that span 
orders #9-14 (^ 4300-7000 A). The spectrum in each order 
was normalized by a fifth-order polynomial. In the analysis 
that follows, we mainly use the data in these orders because of 
their good signal-to-noise {S/N) ratio and general lack of sky 
emission lines and telluric absorption features (relative to red¬ 
der orders). The most prominent features in the spectrum are 
the broad, double-peaked Balmer lines Ha A6563, H)3 /14861 
and H 7 A4341, which are the canonical signatures of a black- 
hole accretion disk. Also present are broad disk lines of He, 
notably He I A5875. Other features such as Na D AA5890- 
5896, diffuse interstellar bands (DIB) at A5780, A6280, and 
telluric features redward of A6850 are also evident. 

2.2. Optical Photometry 

Photometric data were obtained during the same observing 
run using the du Pont 2.5 m telescope, which is also located at 
ECO. Seventy-two V-band images were acquired; by design, 
each 10 -min exposure was strictly simultaneous (time differ¬ 
ence < 1 s) with its corresponding MagE spectroscopic obser- 
vatioiQ. The typical seeing was Ri 0.7" (with a range of 0.6- 
1.1"). We used the 2048 x 2048 Tek#5 CCD, which provided 
a field of view of 8.85' x 8.85' at 0.259" per pixel. Bias cor¬ 
rection and flat-fielding were performed with standard IRAF 

See http://code.obs.camegiescience.edu/mage-pipeline 

The only exception is the first image whose exposure time was only 60 
seconds. 
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Fig. 1.— The mean Magellan/MagE spectram of NovaMus in the observer’s rest frame, covering MagE orders #9-14 (4000-7000 A). Broad and double-peaked 
Balmer lines (Hce, H/1, Hy, and H5) are prominent in the spectrum. He I A5875 emission is also evident. The weak, broad emission features redward of Hyand 
Ha are likely to be He I A4471 and A6678, respectively. The Na D AA5890-5896 doublet, DIB features at A5780 and A6280 (which also has a strong telluric 
component), as well as the telluric feature at A6860 are also labeled. The NovaMus spectrum covering the full MagE wavelength range (~ 3000-10000 A) is 
available in the online journal in electronic foiTn. 


tasks (zerocombine, flatcombine, and ccdproc). The 
centroid of the optical counterpart was determined using the 
CIAO routine wavdetect (Freeman et al. 2002), and the 
quality of the centroiding was confirmed visually. Aperture 
photometry was performed using the IDL procedure aper. 
The radius of the aperture was fixed to eight pixels (2.1"), 
which is the largest aperture that securely excludes significant 
contamination of NovaMus by a faint field star. Choosing a 
smaller aperture (e.g., ~ 5 pixels) would provide very similar 
results. We first produced a light curve of NovaMus relative 
to nearby field stars and then derived an absolute calibration 
relative to a pair of nearby (24" and 41") standard stars. Our 


estimate of the uncertainty in the zero-point is 0.05 mag. 


3. RADIAL VELOCITY ANALYSIS 

The radial velocity K 2 of the secondary star was determined 
by sequentially cross-correlating the 72 object spectra with 
the spectrum of a template star. This was done in turn for in¬ 
dividual echellette orders, rather than using a single merged 
spectrum. The best quality results were obtained for order 
#12 (covering ^ 4700-5500 A), as expected since strong stel¬ 
lar absorption features are present in this order (e.g., the Mg 
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TABLE 1 

Radial Velocity Analysis Results 


Order # 

K2 

(kms“‘) 

7 

(km s“’) 

7b-2454900 
(d) 

No. of 

Object Spectra 

zVv 

A Coverage 

(A) 

Masked Features 

9 a 

413.6±4.3 

17.4±3.1 

46.90550 ± 0.00072 

39 

1.55 

6300-7300 

Ha A6563, 6820-7020 A, 7150-7200 A 

10 

409.5 ±2.7 

19.4±2.0 

46.90351 ±0.00044 

72 

0.92 

5680-6630 

NaD AA5890,5896, DIB A6280 

01 A6300, Ha A6563 

11 

413.0±3.1 

15.2±2.3 

46.90258 ±0.00051 

67 

0.74 

5170-6000 

0 I A5577, Na D AA5890,5896 

12 

401.2 ±2.6 

9.7±1.9 

46.90327 ±0.00044 

72 

0.87 

4700-5500 

H/3 A4861 

13 

407.0±4.2 

17.0±3.1 

46.90211 ±0.00069 

57 

0.92 

4370-5100 

up A4861 

14 

403.9 ±3.9 

9.2 ±2.9 

46.90086 ±0.00068 

63 

1.02 

4060^730 

H 7 A 4341 

Average 

406.8 ±2.2 

14.2±2.1 

46.90278 ±0.00042 






Note. — The quoted uncertainties ai'e at the 1 (7 level of confidence. 

^ The best-fit parameters for order #9 were not included in calculating the average parameter values; for this order, the wavelength ranges 6820- 
7020 A and 7150-7200 A were masked to avoid strong telluric absorption features and sky emission lines. 
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Fig. 2.— The average TDR value of the cross-con'elation between the spec¬ 
trum of each template stai* and NovaMus for order #12, plotted against the 
B — V color of the star. The template star with the highest average TDR value 
is HD 170493, with B — V — 1.11, which coiTesponds to a spectral type of 
r^K5 with Teff ~ 4400 K, based on the scale of S — V color vs. Tg^f vs. Spec¬ 
tral type, which is shown at the top of the figure. 


features at '^5200 A region), and it is relatively free of con¬ 
taminating telluric and interstellar features, 

3,1, Choosing the Best Template Spectrum 

The cross-correlation analysis was performed using the 
software package fxcor as implemented in the IRAF pack¬ 
age, The task fxcor returns the parameter R defined by Tonry 
& Davis (1979; TDR hereafter) provides an estimate of S/N: 
the higher the TDR value, the better the quality of the cross¬ 
correlation, From among our dozens of template spectra, we 
chose the one that gives the highest TDR value, while focus¬ 
ing on the NovaMus data for order #12, Specifically, we: (1) 
selected a K-type template spectrum that gives relatively high 
values of TDR and cross-correlated it against the 72 spectra of 
NovaMus; (2) we velocity-shifted all the NovaMus spectra to 
the rest frame of the template star and summed them; (3) we 
then cross-correlated this summed spectrum against all of our 
template spectra in turn and selected the template spectrum 
that yielded the maximum TDR value. 

The results of this analysis are summarized in Fig,|2l which 
shows how TDR varies among the template stars, which we 
have ordered by their B — V color given in the SIMBAD 
database; we regard the photometric color as a more accurate 
measure of the effective temperature than the spectral type 


(that is given in the same databaseQ The star at the peak 
of this ordered distribution and having the maximum TDR is 
FID 170493 with B — V = 1,11, corresponds to an effective 
temperature Teff ~ 4400 K {Allen’s Astrophysteal Quantities, 

4th edn.). We therefore choose HD 170493 as the master tem¬ 
plate for our cross-correlation analysis. Potential syste matic 
effects related to this choice of template are discussed in il 6 ,ll 

3,2, Fitting the Radial Velocity Curve 

With HD 170493 as the template, we first cross-correlated 
each of the 72 spectra for order #12, For both the object and 
template spectra, we masked out the region around the Hj3 
emission line, and also ~ 100 A at both ends of each order 
because the S/N in these regions is poor. We furthermore re¬ 
quired that the TDR value for a correlation exceed 2,5 (similar 
to the criterion in Orosz et al, 1996), All 72 spectra for or¬ 
der #12 delivered reliable measurements of radial velocity, as 
illustrated in the middle-right panel of Fig, [3 We fitted the 
radial velocity curve with the following model, 

F(f) = 7+^2Cos(27r (2) 

where 7 is the systemic velocity of in the heliocentric frame, t 
is the observation time in Heliocentric Julian Days (HJDs), 
and To is the time of maximum radial velocity; the ra¬ 
dial velocity we used for the template star HD 170493 is 
7 o = —55,07 ±0,07 km s^5 We initially fixed the orbital pe¬ 
riod to the spectroscopic value given by Orosz et al, (1996), 

P = 0,4326058(31) d. Then, using our value of Tq and the 
one in Orosz et al, we obtained a refined value of the orbital 
period: P = 0,43260249(9) d. This period and other best-fit 
parameters are listed in Table [T] 

As a check on our value for the Jif-velocity, we performed an 
independent analysis using the xcor procedure in the MOLLY 
package developed by T, MarshHin place of fxcor in IRAF, 

We again analyzed only the data for order #12 and used the 
same template spectrum (HD 170493) as before. Both the 
object and template spectra were rebinned to the same veloc¬ 
ity dispersion using the vbin procedure, and were normal¬ 
ized using a fifth-order polynomial, as required by the xcor 
procedure. Fitting a Gaussian to the profile of the cross¬ 
correlation function using the MOLLY procedure mgf it, we 

We assume that the reddening of these template star can he ignored due 
to their proximity to the Sun (< 100 pc). 

See http://deneh.astro.warwick.ac.uk/phsaap/software/molly/html/INDEX.html 
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Fig. 3.— The radial velocities of NovaMus relative to the spectral template HD 170493, in MagE orders #9-14. The best-fit radial velocity curves and residuals 
for each order are shown. 
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TABLE 2 

Orbital Parameters for NovaMus 


Parameter 

Value 

Orbital period P (days) 

0.43260249 ± 0.00000009 

K 2 velocity (km s“^) 

406.8 ±2.7 

7 velocity (km s“ ^) 

I4.2±6.3 

To, spectroscopic (HJD—2454900) 

46.90278 ±0.00062 

Mass function /(M) (Mq) 

3.02 ±0.06 

vsini (km s“^) 

85.0±2.6 

Mass ratio q 

0.079 ±0.007 


Note. — The quoted uncertainties are at the IfJ lev el of confi¬ 
dence and include an allowance for systematic en'or; see 


obtained values for both the /T-velocity and its uncertainty 
{K 2 = 402.9 ±1.3 km s^*) that are consistent with those ob¬ 
tained using the IRAF routines. 

In order to obtain a more precise value of K 2 , we analyzed 
data for orders #10 and #11 (which lie to the red of order #12) 
and orders #13 and #14, again masking out ^ 100 A at both 
ends of each order; the wavelength ranges of the orders are 
given in Table [T] As before, we masked out disk emission 
lines (mainly Balmer lines), and we also masked out a num¬ 
ber of strong telluric and interstellar features, which are listed 
in Table [D. As was the case for order #12 discussed above, 
all 72 object spectra for order #10 yielded reliable measure¬ 
ments of radial velocity (i.e., TDR > 2.5); the radial velocity 
curve for this order is shown in Fig. |3] For the other three 
orders, some object spectra were rejected because: (1) the 
TDR value for the correlation was <2.5; (2) the fit to the 
cross-correlation peak did not converge; or (3) the value ob¬ 
tained for the radial velocity lies more than >3(7 away from 
the fitted curved The velocity data for each order individu¬ 
ally were fitted to the same model used for order #12 (Eqn. 2). 
For each of the five orders. Table [T] lists the best-fit parame¬ 
ters, values of reduced and the number of useful spectra. 
Using precisely the velocity uncertainties returned by f xcor, 
we obtain good fits for all five orders, with values that 
are close to unity. We note that there are small differences 
among the best-fit parameters at the level of ^ 'i-Aa level. 

Table |2] provides a summary of the orbital parameters. 
The error-weighted mean of K 2 among the five orders is 
406.8 ±2.2 km s~ ^ A fter including our estimate of the sys¬ 
tematic error (see 1 16.1b . we obtain for the radial velocity semi¬ 
amplitude K 2 = 406.8 ±2.7 km s^*. Using this value of K 2 , 
the orbital period P (Table |2]i and Eqn. (1), we derive a pre¬ 
cise value for the mass function of f{M) = 3.02 ±0.06 Mq, 
thereby confirming the black hole nature of the compact ob¬ 
ject. 

In addition to orders #10-14, and as discussed in the fol¬ 
lowing section, we also performed a radial velocity analysis 
of order #9 for the purpose of measuring the mass ratio q and 
investigating the effects of disk veiling. The best-fit parame¬ 
ters for order #9 are included in Tablell] However, we do not 
include the results for this order in our final determination of 
the radial velocity parameters because relatively few spectra 
generated useful velocity data (39 out of 72) and because the 
quality of the fit is relatively poor (X^/v = 1.55). 


Only five data points in all five orders lie >3—4fT away from the fitted 
curve; removing these radial velocity points does not significantly change the 
fitted values of the parameters. 


4. ROTATIONAF VEFOCITY AND DISK VEIFING 
MEASUREMENTS 

4.1. Rotational Broadening and Mass Ratio 

As the next step toward measuring the mass of the black 
hole, we obtain a precise estimate of the mass ratio q. For 
a short-period, mass-exchange binary like NovaMus, one can 
quite confidently assume that the system is tidally-locked and 
that the secondary fills its Roche lobe (Wade & Horne 1998). 
In this case, q can be simply determined by measuring the ro¬ 
tational broadening of the stellar photospheric lines. Specif¬ 
ically, one measures the radial component of the rotational 
velocity of the secondary star, vsin; and uses the following 
equation (Wade & Horne 1998), 

^ = 0.462,>/3(i+,)2/3, (3) 

Ri 

In measuring vsin/ we use the method of optimal subtrac¬ 
tion. Fig. |4] illustrates the principle of the method using the 
data for orders #10 and #12. Before doing the optimal sub¬ 
traction, we first normalized both the object and template 
spectra using a 5th-order polynomial. The spectra of Nova¬ 
Mus contains contributions from both the secondary (numer¬ 
ous rotationally-broadened absorption features) and the accre¬ 
tion disk (broad disk emission lines like Ha). The best-match 
template spectrum will have virtually the same set of absorp¬ 
tion features as NovaMus; however, the features in the object 
spectrum will be shallower because they are diluted by the 
disk continuum emission. In obtaining the best match of a 
given template spectrum to the object spectrum, one repeti¬ 
tively performs a pair of operations. First, one broadens the 
narrow-line template spectrum, and then multiplies the spec¬ 
trum by the factor /star (0 ^ /star ^1), which represents the 
fraction of the light contributed by the secondary. Finally, one 
subtracts the trial template spectrum from the spectrum of No¬ 
vaMus and determines quantitatively (see below) the template 
spectrum that is freest of photospheric absorption features. 
The results of such an analysis for two orders are shown in 
Fig .|4] The residual spectrum for order #10 (upper panel) only 
shows the broad Ha disk line and the DIB at A6280. The only 
clear ^ectral feature in the residual spectrum of order #12 
(Fig. Hi lower panel) is the Fe II multiple! 42 at /4924, /5018, 
and A5169 (Moore 1972), a feature that is also present in the 
spectra of A0620—00 (Marsh et al. 1994) and the neutron-star 
transient EXO 0748—676 (Ratti et al. 2012). 

In practice, we use the spectrum of our velocity-template 
star HD 170493 and the phase-averaged spectra of NovaMus 
for each order to increase the S/N. Before averaging the spec¬ 
tra for a given order, we shifted the individual spectra to the 
rest frame of HD 170493 using the data in Table [T] We first 
broadened the template spectrum using a grid of trial values 
from 5 to 300 km s^^ in steps of 5 km s^^ We also arti¬ 
ficially smeared the template spectrum using the appropriate 
phase and integration time in order to simulate the velocity 
smearing during an observation of NovaMus. A x'^ test was 
performed on the residual spectrum to determine a minimum 
and the corresponding approximate value of vsin i. Now, start¬ 
ing with this approximate value, we repeated the procedure 
with a trial grid extending from 60 to 120 km s^^ in steps of 
1 km s^*. Making use of the rbroad and optsub tasks in the 
MOLLY, curves of vsin/ vs. x'^ were generated (see Fig. |5]l 
and fitted with a 3rd-order polynomial; our adopted value of 
vsin / corresponds to the minimum of the fitted curve. All the 
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X ( A) 


Fig. 4.— The results of applying the optimal subtraction procedure for order #10 (upper panel) and order #12 (lower panel). Each panel contains the average 
normalized spectrum of NovaMus, the normalized spectrum of the template star HD 170493, and the residual spectrum after optimal subtraction. The Ha 
emission line and strong DIB feature are labeled for order #10. The upper panel also contains the averaged Ha emission line in the observer’s frame, which 
shows the characteristic double-peaked profile. The Fe II multiple! 42 features in the residual spectrum of order #12 are labeled in the lower panel. All the spectra 
are smoothed with a 5-pixel boxcar for the puipose of presentation. 
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Fig. 5.— Measurement of the rotational velocity vsini in MagE orders #9-14. For each order, the upper panel shows the curve of vsin; vs. which we 
fitted with a cubic polynomial. Our adopted value of vsini is the value corresponding to the minimum in x^- The lower panel illustrates the method we used to 
estimate the uncertainty in v sin i via the bootstrap method described in the text. The histogram of v sin i is fitted with a Gaussian in determining the uncertainty. 


emission lines, interstellar features, and telluric features are 
masked out during the optimal subtraction. 

The effect of limb darkening needs to be properly consid¬ 
ered in the optimal subtraction. Since this effect varies with 
wavelength, we calculated a linear limb darkening coefficient 
for each order of the MagE spectra based on the data given 
by Claret et al. (2012). We retrieved from their table the limb 
darkening coefficients for the UBVRI bands. Eor the effec¬ 
tive temperature of the companion star, we adopted the value 
Teff = 4400 K. The star’s surface gravity log,g (in cm^^ s^*) 
is securely in the range 4.0^.5, and for each band we adopted 
the average of the tabulated values of the limb darkening coef¬ 
ficient for logg = 4.0 and \ogg = 4.5. Einally we linearly in¬ 
terpolated between the UBVRI bands using the central wave¬ 
length of each order, thereby obtaining the values we adopt 
for the limb darkening coefficients (see the third column in 
Table [2i, which range from 0.71 (Order #9) to 0.91 (Order 
#14). 


We utilized the bootstrap method to determine the uncer¬ 
tainty in vsin/ following Steeghs & Jonker (2007) and Ratti 
et al. (2013). We generated 500 simulated NovaMus spectra 
using the MOLLY task boot. The bootstrap method randomly 
selects data points from the original spectrum while keeping 
the same total number of data points in the spectrum. Eor each 
simulated spectrum, vsin/ was measured following the same 
procedure as described above. The 500 vsin/ values follow 


TABLE 3 

Limb Darkening Coefficient and Rotational Velocity 


Order # 

Central X (A) 

u (limb darkening) 

vsin/(kms 

Disk Veiling (%) 

9 

6800 

0.706 

83.8±7.5 

40.3 ± 3.4 

10 

6130 

0.762 

87.2±3.5 

5L4±L3 

11 

5590 

0.808 

88.6±4.5 

56.7 ±1.4 

12 

5140 

0.845 

82.8 ±3.5 

59.3 ±1.0 

13 

4750 

0.878 

79.4±5.9 

55.8±2.1 

14 

4400 

0.908 

85.5±7.1 

69.1 ±1.5 

Average 



85.0±L3 



Note. — The quoted uncertainties are at the IfT level of confidence. 


a Gaussian distribution (see Eig. |5ll. The mean of the Gaus¬ 
sian distribution is very close (within 1 km s^*) to the fitted 
values of vsin/ discussed above. Eor each order, we take the 
mean of the Gaussian distribution as the final measurement of 
vsin/ and the standard deviation of the mean a as the uncer¬ 
tainty. The measured values of vsin/ for six orders (#9-14) 
are presented in Eig.|5]and TableO. 

The error-weighted mean vsin / of all six orders is (vsin/) = 
85.0 ± 1.3 km s^*. After including an allowance for sys¬ 
tematic uncertainty, related to the choice of te mpla te spec¬ 
trum and the limb darkening coefficient (see 116. Il l, our fi¬ 
nal adopted value is vsin/= 85.0±2.6 km s^^ Using this 
value and our fC-velocity for the secondary {K 2 = 406.8 ± 
2.7 km s^*; 113.21) . and Eqn. (3), we obtain for the mass 
ratio the value q — 0.079 ± 0.007. Using this mass ra- 
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veiling factor 


Fig. 6 .— Measurement of the disk veiling fraction in order #9-14. In each 
panel, the dotted line indicates the disk veiling fraction (labeled by the “mean 
veiling 1” value) of the averaged spectrum. The histogram shows our disk 
veiling measurements for individual spectra; the dashed line indicates their 
mean value as obtained from the histogram (labeled by the “mean veiling 2” 
value). 

tio and the value of the mass function given in il3.21 the 
masses of the black hole and secondary star are M = 3.52 ± 
0.07(sin^^ i) M© and M 2 = 0.28 ± 0.03 (sin^^ i) M©, respec¬ 
tively. 

4.2. Disk Veiling 

The disk veiling factor /disk; i e., the fraction of the total 
light that is non-stellar, is also determined by the optimal sub¬ 
traction method. For each order, we hrst broaden the HD 
170493 te mpla te spectrum using the values of vsini deter¬ 
mined in Mill and then optimally subtract it from the spec¬ 
trum of NovaMus. The procedure gives the fraction of light 
due to the secondary star, /star, for the template spectrum that 
corresponds to the minimum value of x^- The disk veiling 
factor is then simply /disk = 1 — /star- As in determining v'sin/ 
we computed /disk using for each order the averaged, rest- 
frame spectrum of NovaMus. The values we adopt for the 
disk veiling factor for the six orders are indicated in Fig.|6]by 
dotted lines and are listed in the last column of Tabled 

As a check on the robustness of our measured values of 
/disk, we repeated the optimal subtraction procedure for each 
spectrum individually, using only those spectra that produced 
a reliable measurement of radial velocity (see A his¬ 
togram of /disk for each order is shown in Fig. |6l the dashed 
line in each panel represents the weighted mean value of the 



A ((i m) 


Fig. 7.— The simplistic spectral energy distribution of NovaMus, decom¬ 
posed into the secondary star component (red line) and the accretion disk 
component (blue line). The pink line shows the sum of the two components. 
The solid blue and pink lines cover the observed wavelength range for MagE 
orders #9-14. The dashed lines are extrapolations to longer wavelengths. The 
star component is represented by a 4400 K blackbody spectnim. The blue line 
represents a power-law model obtained by fitting the observed relative flux of 
disk emission (filled blue squares) in the six orders, excluding the data point 
for order #13 (the second blue square from the left),, which is a cleai* outlier. 
The black dotted lines show the effective wavelengths of the R, I, J, H, K 
bands. Our estimate of the disk veiling factor in each band is also labeled. 

disk veiling factor. These alternative measurements of /disk, 
are consistent within ~ 1 O’ (for all orders) with the values de¬ 
termined as described above using the phase-averaged spec¬ 
tra. We note that for hve of the six orders /disk is greater than 
0.5, i.e., the disk is brighter than the secondary star. 

Fig. |6] indicates that the disk veiling factor generally de¬ 
creases with increasing wavelength, i.e., the disk veiling is 
less signihcant for redder wavebands. However, this does not 
necessarily mean that disk veiling can be ignored, even in the 
J and K bands, as claimed by Gelino et al. (2001) and oth¬ 
ers. We utilize a simplistic spectral energy distribution (SED) 
model and an extrapolation to illustrate this point (see Fig.lT]). 
Our crude SED by no means represents an accurate model 
of the disk emission. Eor the stellar component, we assumed 
blackbody emission at 4400 K (red line in Eig.|2l see ij3]and 
Fig. El), and we then computed the relative disk flux in each 
order (hlled blue squares) based on our measurements of /disk- 
We htted these estimates of disk flux using a power-law model 
to represent the SED of the disk component and then extrapo¬ 
lated the power law into the near-infrared (dashed blue line). 

At shorter wavelengths, the disk component does indeed 
decrease rapidly with increasing wavelength so that the stel¬ 
lar component becomes dominant as one approaches the R- 
band. However, at longer wavelengths the blackbody com¬ 
ponent begins to fall rapidly, and disk veiling reaches a min¬ 
imum (Ki 40%) at around 1 fim and then begins to increase. 
Eor our data and this crude model, we estimate that /disk for 
the /, J, H, K bands is in the range of 40-60% (see the black 
dotted lines in Eig. El*, which is consistent with the estimates 
of Reynolds et al. (2008) for NovaMus of 40-50% in the near- 
infrared. We note that /disk in the near-infrared could be lower 
than our estimate if the disk component of emission has a 
thermal spectrum that falls below the extrapolated power-law 
spectrum, which we assume. However, we conclude that it is 
not warranted to assume that the effects of veiling are negli¬ 
gible in the near-infrared. 
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Fig. 8 .— The original (left panel) and phase-binned (right panel) V-band light curves of NovaMus, acquired with the du Pont 2.5-m Telescope. In the left panel, 
the open (filled) circles are the photometry data points obtained during the first (second) night of observation. The error bar in the lower left corner represents the 
typical uncertainty. In the right panel, the light curve is binned into 20 phase bins, each of which has three or four individual data points. 



5. LIGHT CURVE 

Our calibrated and phase-folded V-band light curve, which 
was obtained in exposures that obt ained simultaneously with 
the spectroscopic exposures (see 112.2b . is shown in Fig. 0 
Zero photometric phase, which is equivalent to 0.75 spectro¬ 
scopic phase, corresponds to the inferior conjunction of the 
secondary. Two versions of the light curve are shown, un¬ 
binned and binned into 20 phase bins, each with 3-4 data 
points per bin. 

During our observations, NovaMus was ^ 0.5 magnitude 
brighter on average in the V-band than it was in earlier pho¬ 
tometric monitoring campaigns (e.g., during 1992-1995, as 
reported by Orosz et al. 1996). Our V -ban d light curve shows 
substantial aperiodic flickering (see 1 16.2b . Nevertheless, the 
phase-binned light curve exhibits the characteristic signatures 
of ellipsoidal modulation, with two minima at phases ~ 0.0 
and ^ 0.5 and two maxima at phases ^ 0.25 and 0.75. The 
two maxima do not have equal amplitude as one would expect 
for an ideal ellipsoidal modulation, which could be caused by 
a hot spot in the accretion disk. Such unequal maxima have 
frequently been observed in the light curves of other black 
hole binaries in quiescence, e.g., A0620—00 (McClintock & 
Remillard 1986; Cantrell et al. 2010). 

6 . DISCUSSION 

6.1. Assessment of Systematic Uncertainties 

We identify sources of systematic error, and obtain esti¬ 
mates of their magnitude, first for the orbital parameters K 2 
and Tq and the systemic velocity 7 , and then for the rotational 
velocity vsin/. 


6.1.1. Orbital Parameters and Systemic Velocity 

We first consider the uncertainty associated with the choice 
of template spectrum. While we chose the spectrum of 


HD 107493 because it gave the maximum TDR value ( il3.1b . 
obviously this spectrum is not identical to that of the Roche- 
lobe-fllling companion star. Focusing on the high-quality data 
for order # 12 , we repeated our analysis for six other template 
spectra with TDR values > 60 (see Fig. |2l); the results are 
summarized in Table |4] The spread in the values of K 2 is 
< 0.8 km s^*, which is significantly smaller than the statisti¬ 
cal uncertainty. As our estimate of the systematic error in K 2 
associated with the choice of the template, we take the stan¬ 
dard deviation of the seven values of K 2 , which is 0.3 km s^*. 
In like manner, we obtain estimates of systematic error of 
4.2 km s^^ for 7 , and of 0.00013 day for Tq. 

A second source of systematic uncertainty is the effect of 
the different widths of the photospheric lines for NovaMus 
which are rotationally broadened, and for the slowly-rotating 
template star, which are narrow. To assess this effect, we 
broadened the lines of our chosen template (HD 107493) to 
85 km s^^ (1 14.lb and repeated the radial velocity analysis. 
The change in K 2 is 0.5 km s^* (Table|4]i, which we adopt as 
our estimate of systematic error; we obtain similar estimates 
of systematic error for 7 and Tq of 0.8 km s^^ and 0.00020 
day, respectively. 

For each of the three parameters in question, we obtain our 
final estimate of error, which is given in Table 2, by linearly 
adding the larger of the two systemat ic er rors discussed above 
onto the statistical error obtained in 1 13.2l (given in Table 1). 


6.1.2. Rotational Velocity 

We identify two sources of systematic uncertainty for v sin i, 
namely, the value adopted for the limb darkening coefficient 
and the choice of template spectrum. Again, using the data 
for order # 12 , we considered an uncertainty of ± 0.1 in our 
adopted value of 0.845 (Table 3) and And that this corresponds 
to an uncertainty of < 1.5 km s^^ (see Table |5]). We then 
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TABLE 4 

Radial Velocity Analysis Results eor Different Templates (Order #12) 


Templates 

B-V 

Template 7 

TDR 

Ki 

7 

To - 2454900 

No. of 

zVv 

(Order #12) 


(km s“’) 


(km s“’) 

(km s *) 

(d) 

Object Spectra 


HD 147776 

0.95 

7.4 ±0.1 

60.8 

402.0 ±2.7 

17.1±2.0 

46.90339 ±0.00044 

72 

0.83 

HD 130992 

1.02 

-57.1 ±0.1 

64.1 

40L9±2.7 

10.9 ±2.0 

46.90333 ±0.00044 

72 

0.79 

HD 31560 

1.09 

6.4 ±0.1 

74.3 

40L6±2.7 

11.4±2.0 

46.90332 ±0.00044 

72 

0.81 

HD 170493 

1.11 

-55.1 ±0.1 

78.4 

40L2±2.6 

9.7 ±1.9 

46.90327 ±0.00044 

72 

0.87 

HD 170493 (broadened) 

1.11 

-55.1±0.1 

76.7 

40L7±3.0 

10.5 ±2.2 

46.90347 ±0.00050 

72 

0.60 

HD 131977 

1.11 

26.8 ± 0.1 

65.9 

401.6 ±2.7 

10.9 ±2.0 

46.90326 ±0.00045 

72 

0.79 

HD 156026 

1.16 

0 . 1 ± 0.1 

61.6 

40L9±2.8 

4.6 ±2.0 

46.90344 ±0.00046 

72 

0.87 

HD 120467 

1.26 

-37.8±0.1 

63.1 

40L9±2.8 

5.5±2.1 

46.90365 ±0.00047 

72 

0.79 


TABLES 

Rotational Velocity for Different Templates 
(Order #12) 


Template 

B-V 

LI (limb darkening) 

vsin/(kms 

HD 170493 

1.11 

0.845 

82.8±3.5 

HD 170493 

1.11 

(1745 

8I.7±3.6 

HD 170493 

1.11 

0.945 

84.3 ±3.7 

HD 147776 

0.95 

(1^23 

82.9±4.6 

HD 130992 

1.02 

0.842 

82.3 ±4.5 

HD 31560 

1.09 

0.845 

82.4±4.2 

HD 131977 

1.11 

0.845 

82.3 ±4.0 

HD 156026 

1.16 

0.835 

83.8±4.3 

HD 120467 

1.26 

0.796 

85.5±5.2 


Note. — The quoted uncertainties are at the Iff level of 
confidence. 


compared values of vsini derived using once again the six 
template spectra with TDR values > 60. For each template 
spectrum, we used the appropriate limb darkening coefficient 
(Claret et al. 2012) based on the B — V color and correspond¬ 
ing effective temperature of the star. The standard deviation 
for the six values of vsin; is 1.2 km s^^ Using these two esti¬ 
mates of systematic error as a guide, we allow for systematic 
error by simply doubling the statistical error bar obtained in 
il4. 11 (1.3 km s^*). Our final adopted value of vsin/ given in 
Table 2, is therefore 85.0±2.6 km s^* 

6.2. Comparison to Previous Work 

We have obtained precise measurements of the mass func¬ 
tion f{M) and mass ratio q for NovaMus using multi-order 
echellette spectral data obtained with Magellan/MagE. Com¬ 
pared to the data analyzed in previous work (e.g., Orosz et al. 
1996 and Casares et al. 1997), our spectra have higher spec¬ 
tral resolution (by a factor of ^ 3-5) and provide more com¬ 
plete phase coverage. We also observed many more template 
stars (38), all of them during the same observing run, which 
has enabled us to choose a template spectrum that is an excep¬ 
tionally close match to the spectrum of NovaMus. Analyzing 
data independently for several orders is also a feature of our 
work. All of these beneficial factors have contributed to qual¬ 
ity of our work. 

Our value for the radial velocity amplitude {K 2 = 406.8 ± 
2.7 km s^*) is very close to that obtained by Orosz et al. 
(1996; K 2 = 406 ±7 km s^*), while the uncertainty has been 
reduced significantly. Meanwhile, our K 2 value is barely 
consistent with that of Casares et al. (1997; K 2 = 420.8 ± 


6.3 km s^'), differing by Ri 2a. We have substantially in¬ 
creased the precision of the orbital period (by a factor of > 30) 
by extending the baseline of observations from those made 
earlier (Remillard et al. 1992; Orosz et al. 1996) out to our 
observations in 2009. 

Our mass ratio {q = 0.019 ± 0.007) agrees with the approx¬ 
imate value (0.0476 < q < 0.0833) found by Antokhina & 
Cherepashchuk (1993), which was obtained by modeling an 
/-band light curve. It is significantly smaller than the less re¬ 
liable estimate of ^ = 0.133 ±0.019 obtained by Orosz et al. 
(1996; see ©■ Casares et al. (1997), who used the same 
direct method of measurement we used, likewise obtained a 
mass ratio greater than ours, o' = 0.13 ± 0.04. (Their estimate 
of the rotational broadening isvsin/=106±13 km s ^ ^ com¬ 
pared to our value of vsini = 85.0±2.6 km s^*.) The relative 
virtues of our measurement is our better spectral resolution 
(« 60 km s^^ at ~ 5000 A) and our ample selection of tem¬ 
plate spectra, which allowed us to closely match the spectrum 
of the donor star; Casares et al. had applied a single template 
spectrum of a KOV star. 

An important feature of our work is the simultaneity of 
the spectroscopic and photometric observations, which allows 
one to accurately correct for the effects of disk veiling, de¬ 
spite the aperiodic variability of the non-stellar component of 
light. Unfortunately, during our observation NovaMus was in 
an “active” quiescent state (Cantrell et al. 2008), and the light 
curve is strongly affected by aperiodic flickering and does not 
provide a useful constraint on the inclination angle. Neverthe¬ 
less, our simultaneous observations provide very secure mea¬ 
surements of the disk veiling factor at six wavelengths. In a 
subsequent paper, we will use these data in conjunction with 
light curves available in the literature, which were obtained in 
the “passive” quiescent state (Cantrell et al. 2008) to constrain 
the inclination angle and black hole mass. 

6.3. The Li X6708 Feature 

The Li /16708 absorption line was previously detected in the 
spectra of NovaMus (Martin et al. 1996). This feature also ap¬ 
pears in the spectra of other similar black hole soft X-ray tran¬ 
sients like A0620-00 (Marsh et al. 1994) and GS 2000±25 
(Filippenko et al. 1995). This feature is potentially relevant in 
exploring the Li production mechanism around a black hole 
(Martin et al. 1994). Meanwhile, the variation of the line in¬ 
tensity with orbital phase may reveal effects of irradiation on 
the formation of the Li feature (Martin et al. 1996). 

We find that the Li /16708 line is marginally detected in the 
phase-averaged spectrum of NovaMus and also in the residual 














12 


WU ET AL. 



X ( A) 


Fig. 9.— The averaged spectrum of NovaMus (top), the template spectrum of HD 170493 (middle), and the residual spectrum (bottom) after optimal subtraction, 
all for MagE order #9 and in the heliocentric frame. The location of Li A6708 is indicated by the dotted line. This Li absorption feature is marginally detected 
(~ 3.5<J). The inset plot shows the phase-resolved NovaMus spectra around Li A6708. The Li feature is detected in the lower spectrum, in which the photometric 
phase 0 is around zero (0 > 0.75 or 0 < 0.25). There is no cleai* Li feature in the upper spectrum ai'ound photometric phase 0.5 (0.25 < 0 < 0.75). All the 
spectra are smoothed with a 5-pixel boxcar for the purpose of presentation. 


spectrum after optimal subtraction (Fig.|9]l. We find an equiv¬ 
alent width of EW(Li)= 210 ±60 mA, after correcting for the 
disk veiling in order #9, for the NovaMus spectrum averaged 
over the full orbital phase. Following Martin et al. (1996), we 
investigated the variation of the EW(Li) with orbital phase 
(see the inset plot of Fig. IS. At zero photometric phase 
(0 > 0.75 or 0 < 0.25), we obtained EW(Fi)= 550 ±90 mA, 
which is slightly higher than the value in Martin et al. (1996; 
420 ± 60 mA). The Fi A6708 feature is not detected at pho¬ 
tometric phase 0.5 (0.25 < 0 < 0.75), which agrees with the 
findings in Martin et al. (1996). We note that our EW(Fi) 
for NovaMus is slightly higher (within la) than the value re¬ 
ported for A0620—00 of EW= 160 ± 30 mA by Marsh et al. 
(1994). 
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